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Abstract. We construct anisotropic conductivities with the same Dirichlet- 
to-Neumann map as a homogeneous isotropic conductivity. These conduc- 
tivities are singular close to a surface inside the body. 



1 Introduction 

An anisotropic conductivity on a domain 0, C R" is defined by a symmetric, 
positive semi-definite matrix- valued function, a = {a'^^{x)). In the absence 
of sources or sinks, an electrical potential u satisfies 

{V-(7V)u = dj(7^''{x)dku = OinO, (1) 
u\dn = f, 

where / is the prescribed voltage on the boundary. Above and hereafter 
we use the Einstein summation convention where there is no danger of con- 
fusion. The resulting voltage-to-current (or Dirichlet-to-Neumann) map is 
then defined by 

A,(/) = Bulan (2) 

where 

Bu = v^a^^dkU, (3) 

u is the solution of (1) and v = (z/i, . . . , z/„) is the unit normal vector of dQ. 
Applying the divergence theorem, we have 
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where u solves (1) and dS denotes surface measure on 90. QaU) represents 

the power needed to maintain the potential / on dVL. By (4), knowing 

Qa is equivalent with knowing Ag-. li F : ^ Vt, F = 

is a diffeomorphism with F\qq^ = Identity, then by making the change of 

variables y = F{x) and setting u = v o F^^ in the first integral in (4), we 

obtain 

^F,a = Act, 

where 

1 ^ f) F^ f) F^ 

is the push- forward of the conductivity a hy F. Thus, there is a large 
(infinite-dimensional) class of conductivities which give rise to the same 
electrical measurements at the boundary. The uniqueness question, first 
proposed by Calderon [C], we wish to address is whether two conductivi- 
ties with the same Dirichlet-to-Neumann map must be such pushforwards 
of each other. By a direct construction, we will show that the answer is no. 
Let D CC O be a smooth subdomain. We will construct in Q a conductiv- 
ity a for which the boundary measurements coincide with those made for 
homogeneous conductivity 7 = 1 in 0. We note that the conductivity is 
singular in the sense that some components of the conductivity tensor go to 
zero, i.e., correspond to perfectly insulating directions and some components 
can go to infinity, i.e., correspond to perfectly conducting directions, as one 
approaches dD. The physical meaning of these counterexamples is evident: 
For instance in medical imaging, e.g., in Electrical Impedance Tomography 
(EIT) certain anisotropic structures can form barriers which, including their 
interior, appear in measurements to be a homogeneous medium. 

To construct these counterexamples, we need to consider a variant of (1), 
which comes from the Laplace-Beltrami operator on a compact Riemannian 
manifold with boundary. Let us assume now that (M, g) is an n-dimensional 
Riemannian manifold with smooth boundary dM. The metric g is assumed 
to be symmetric and positive definite. The invariant object analogous to the 
conductivity equation (1) is the Laplace-Beltrami operator, which is given 
by 

^gU = G-^'^d,{G^l^gi''dku) (6) 
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where G = de\.{gjk), [gjk] = W^] ^- The Dirichlet-to-Neumann map is 
defined by solving the Dirichlet problem 

AgU = in M, (7) 
u\aM = /• 
The operator analogous to is then 

with u = (z^i, . . . , Un) the outward unit normal to dM. In dimension three 
or higher, the conductivity matrix and the Riemannian metric are related 
by 

a^'^ = det((7)^/V^ or ^^''^ = det(a)2/("-2)a^^ (9) 

Moreover, = Ao-, and A^.^ = A^, where ilj*g denotes the pullback of the 
metric (7 by a diffeomorphism of M fixing dM |LeUj . 

In dimension two, (9) is not valid; in this case, the conductivity equation 
can be reformulated as 

Divg(^Gradgn) = in M, (10) 
u\dh[ = f 

where /3 is the scalar function /3 = jdeto"!^/^, g = {gjk) is equal to {(Tjk), 
and Divg and Gradg are the divergence and gradient operators with re- 
spect to the Riemannian metric g. Thus we see that, in two dimensions, 
Laplace-Beltrami operators correspond only to those conductivity equations 
for which det((T) = 1. 

For domains in two dimensions, Sylvester |Sy| showed, using isothermal 
coordinates, that one can reduce the anisotropic problem to the isotropic 
one; combining this with the isotropic result of Nachman |Naj . one obtains 

Theorem 1 If a and a are two -smooth anisotropic conductivities in 
C for which Aq- = A^, then there is a diffeomorphism F : — > 0, 
F'ldn = Id such that a = F^a. 

In dimensions three and higher, the following result is known (see |LU) . 
pTTTT] . and [Lillj): 

Theorem 2 If n > 3 and {M, dM) is a manifold with connected, 
boundary, and g,g are metrics on M such that Ag = Ag, then there 
exists a diffeomorphism F : M ^ M such that F\qd = Id. 
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2 Counterexamples 



Returning now to domains 0, CC R", n > 3, let D CC be an open subset 
with smooth boundary and g = gij be a metric on i}. Let y & Dhe such that 
there is a diffeomorphism F : Q\ {y} ^ Q\D, and let g = F^g onQ\D. To 
obtain a conductivity on all of 0,, first extend ^ to a bounded metric inside 
D and denote this new metric on by g. We make this continuation so that 
the conductivity jumps on dD and that ^ > c > is smooth inside D, but 
is otherwise arbitrary. Let a be the conductivity corresponding to g by (9) 
in We say that v is a solution of the conductivity equation if 

\7-a'Vv{x) = in the sense of distributions in 17, (11) 

v\dn = fo, 

V G L^{n) 

where fo € H2(^dQ). That the equation in the sense of distributions means 
that V G H^{n) and a\/v G H{n;\/-) = {w G L'^{n) : V-w £ ^^(Q)}. 
If this problem has a unique solution, we define the Dirichlet-to-Neumann 
map 

A5/0 = n-aVv\Qn. 
Our aim is to prove the following result, first announced in |GLU| : 

Theorem 3 Let C R", n > 3, and g = gij a metric on Vl. Let D dVi he 
such there is a C°° -diffeomorphism F : Q\{y} — > r2\L> satisfying F\qq = Ld 
and that 

dF{x) > Co/, det{dF{x)) > ci dist^„ (12) 

where dF is the Jacobian matrix in Euclidean coordinates o/R" and cq, ci > 
0. Let g = F^g and g be an extension of g into D such that it is positive 
definite in D*"*. Finally, let'j and a be the conductivities corresponding to g 
andg by (9). Then the boundary value problem for the conductivity equation 
with conductivity a is uniquely solvable and 

As = A^,. 

Note that here is no diffeomorphism H : 0, ^ Q such that a = H^^'j, so the 
Riemannian manifolds corresponding to a and 7 cannot be the same. Also, 
a can be changed in arbitrary way inside D without changing boundary 
measurements. The proof of Theorem IHl will be given below. 
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Example. Let Q = -6(0,2) C be the ball with center and radius 2. 
Consider y = e D = B{0, 1) and the map F : ^} \ {0} ^ Q\D given by 

F:rr^(M + i)^. (13) 
2 \x\ 

Let 7 = 1 be the homogeneous conductivity in and define a = F^^. 
Now the metric tensor g and the corresponding conductivity (Tg are related 
by = g\^/'^ g^^ . Let g be the metric corresponding to 7 and Tj be 
the metric corresponding to cr. Consider these in the standard spherical 
coordinates on 17 \ {0}, {r,(f),9) 1— s- (r sin cos (/>, r sin ^ sin 0, r cos 0) G R^. 
With respect to these coordinates, we see that the metric g and conductivity 
7 correspond to the matrices 

/ 1 \ / ''^sin6l 

5 = , 7 = sin 6* 

\ r2sin2 6' / \ {smO)-^ 

and g and a correspond in the annulus 1 < r < 2 to the matrices 
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Let a be a continuation of a that is C°°-smooth in D. Then these metric 
are as in for Theorem |31 and in particular F satisfies (|12j) . 

To start to prove Theorem |21 we start with the following result: 

Proposition 1 Let Q C R", n > 3, and g = gij a metric on il. Let u 
satisfy 

Agu{x) =0 in Q, 
u\an = /o G C°^{dn). 

Let D d Vl he such that there is a diffeomorphism F : 0,\ {y} Q\ D 
satisfying F\qq = Ld. Let g = F^,g and v he a function satisfying 

Agt;(x)=0 inn\D, 

u\dQ = /o, 

u G L'^inXD). 
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Then u and F*v coincide and have the same Cauchy data on dil., 



duU\dM = (kF*v\QM 



(14) 



where v is unit normal vector in metric g and v is unit normal vector in 
metric g. Moreover, for the constant cq := u{y) we have 



Proof. Let g = {gij) be a Riemannian metric tensor defined on 17 C R". 

First, we continue the C°°-metric gioa, metric 5 in R" such that Qij{x) = 
gij{x) for x G is such that for any y G R™ there is a positive Green's 
function G(x, y) satisfying 



There are several easy ways to obtain this continuation. For instance, we 
can continue g to a metric g such that outside some bah i?(0, R) the metric 
is hyperbohc. This impHes that the manifold (R",g) is non-parabolic and 
has a positive non-constant super-harmonic function (see |Grj ) . By |LiTj 
there thus exists a positive Green's function G{x,y). 

Next we consider the probability that Brownian motion Bf on manifold 
(R",g) sent from the point x at time t = enters an open set U: Let 



Let A = Tl^ \ U . By Hunt's theorem (see, |Hj, or |;Qr|, Prop. 4.4) we have 



where sa{x) is the super-harmonic potential of A, that is, s^(x) is infimum 
of all bounded super-harmonic functions h in R" such that h\u = 1 and 
/i > in R". 

Let U = B{y,r) and 



lim v{x) = Cq. 



A,G{-,y) = 6, 



y 



inR". 



eu{x) = P({there is t > such that Bf G U}) 



eu{x) = sa{x) 



m{y, r) 



inf G{x,y). 

xS:dB(y,r) 



The function 




1 



min {G{x,y),m{y,r)) 



m{y, r) 
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is positive super-harmonic function and satisfies hy^r\B{y,r) = 1- Since 
eu{x) < hy^r{x) and liuir^o m{y , r) = oo, we see tliat 

lim P({there is t > sucli that Bf G B{y, r)}) = 0. (15) 

r— >0 

In particular, taking limit r — > we see that the probability that Bf = y 
for some t > is zero. 

Now, let u be any solution of 

Agu{x) = in 17 \ {y} (16) 

u\dn = fo, 
5GL-(f]\M). 

Denote fr = Sjg^^j^ ,,) and let T{r,x) € (0, cxo] be the first time when Bf G 
9$7 U dB{y,r). Similarly, let t(0, x) be the hitting time to dil.. 

Let next Xo,r be a random variable defined by Xo,r = 1 if Br(r,x) 
and zero otherwise, and let Xir = ^ ~ Xor- Then by Kakutani's formula. 

Letting r ^ and using the fact that | |/r i |oo < 1 1^1 Ico are uniformly bounded 
and ((T3|) we see u{x) = E(/o(i?^j-Q ^p). Thus u{x) = u{x) for x G \ {y} 
where 

Agu{x) = in 17, (17) 
u\dn = fo- 

Thus we have shown that boundary value problem (|l(i|) is solvable, and 
that the solution is unique. Next we change this problem to an equivalent 
problem. Let F : Q \ {y} Q \ D he a diffeomorphism that is identity at 
boundary dQ. Define the metric g = F^,g\ then the boundary value problem 
(|17j) is equivalent to 

Agt;(x) = inn\;D (18) 
v\dQ. = fo, 

V G l°°(o\:d). 

and solutions of the problems (|T7j) and (fT8|) are related by v{x) = u{F{x)) 
for X G r2 \ {y}. 

Clearly, the Cauchy data of the equations (|17j) and (|18j) coincide in the 
sense of ()14[). Moreover, 

lim v{x) = Co := u{y) 

x^dD 
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This concludes the proof of Prop. 1. □ 

We remark that the apphcation of Brownian motion above is not essential 
but makes the proof perhaps more intuitive. Alternatively, in the proof of 
Prop. 1 one can use properties of L^'-Sobolev spaces. Indeed, for u satisfying 
Agu{x) = in Q \ {y}, u\gfi = /o, and u € -L°°(0 \ {y}), we can consider 
the extension u £ L°°(fl) that satisfies 

Agu{x) = F in 17, 

where F is a distribution is supported in y. Then F has to be a finite 
sum of derivatives of the Dirac delta distribution supported at y. Now, 
F G W-^'P{n) for all 1 < p < cx), and since 6y ^ W-^'P{n) for p > ^ 
and the same is true for the derivatives of the delta distribution, we see that 
that F = 0. This implies that u = u. 

We now turn to the proof of Theorem 3. From now on, we assume that 
F -.nXiy} ^n\D, F{x) = {F^{x),..., ^"(x)) is such that condition ^ 
is satisfied. First, continue g into D so that that it is positive definite in 
Next, extend v inside D to a function 

h{x) = v{x) for X G 0, \ D, 
h{x) = Co for x G D. 

Our aim is to show that h is the solution of Since any solution of (|18|) 
is constant on the boundary of dD and g is positive definite in D, we see 
that the solution has to be constant inside D. Thus h is the unique solution 
if it is a solution. 

Now we are ready to show that w is a solution also in the sense of distri- 
butions. First we note that when y{x) = F^^{x) we have v{x) = u{F~^ (x)) 
and 

and since u G H^{n) and G L'^{n\D) we have v G H'^{n\D). 
Also, as h £ C(r2 \ and the trace v ^ v\q£) is a continuous map 

C{n \ Z)™*) nH^{n\D) ^ L'^idD), we see that trace v\3d is weU defined 
and is the constant function having value cq. Since restrictions of /i to D 
and i}\D are in H^{D) and H^{Q\D), correspondingly, and the trace from 
both sides of dD coincide, we see that h G H^{^). 
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Next we show that aVv £ L^{n\D). Let ei, 62, ■ ■ ■ , Gn be the standard 
Euchdean coordinate vectors in R". Then for x € Q\D 

W^'dkvix) = {ej,Vgv{x)yg = {d{F-^)e,ygu{F~\x)))g 

and we see from H12() that this inner product is uniformly bounded. By (jSJ 
and (fT^ . |det(5f)(a;)| < C3 distRn {F^'^{x), y)"^ and thus in Q.\D the functions 
Vj{x) = \deig{x)\^/'^5jkg^^diV satisfy 

\Vj{x)\ < C4 distRn {F-\x),y). 

This imphes that V = {Vi{x), . . . ,Vn{x)) is in H{n\D;V-) n C{Tl\ D'""^). 
Moreover, by !Gr the normal trace W n- W\gD is a continuous map 
H{n\D; V- )n(7(n\L>*"*) ^ H-^/^{dD), we see that n- V\od = 0. Since the 
normal traces of V coincide from both sides of dD we see that V G H{n; V- ) 
and 

n 

^ dk{\detg{x)\^/^g'''dih) =V-V = in Q 
j,k=i 

in the sense of distributions. This means that h satisfies the conductivity 
equation in Q in the sense of distributions. This proves Theorem |3J 

A similar approach can be used to construct a different type of counterex- 
ample. In contrast to the previous one, here the conductivity is not bounded 
above near the singular surface. Consider again the sets O = 5(0,2) and 
D = -6(0, 1) in R^. In spherical coordinates we define the metric g and the 
conductivity a in the domain {1 < r < 2} by the matrices 

/ (r-l)-2 \ 

\ p'^sin^e J 

/ (r - l)p^sme 

a = i (r - l)-isin6l 

\ {r-l)-^sm-^9 

where /) > is a constant. Thus, (fi \ D,g) is isometric to the product 
R-i- X Sp with the standard metric where is a 2-sphere with radius p. 
Now extend the conductivity a to 17 so that it is positive definite in D. It 
can be shown that, in the domain ^}\D, the equation 

S/-aVv{x) = Omn\D 
v\dn = /o, 



9 



has a unique solution. Also, when p is small enough, we can extend the 
definition of v{x) to the whole domain Q by defining it to have the constant 
value Co = limx^dD v{x) in D. The function v{x) thus obtained is a solution 
of the boundary value problem 

V-dVv{x) = in J^, 

v\dn = fo, 

V G L'^in) 

in the sense of distributions. 

In this case, the boundary measurements do not give information about 
the metric inside D. To see this, consider D,g) as R+ x with 
coordinates (t, (p, 9), where t = t{r) = — log(r — 1), i.e., t € R+ and 0, 9 are 
coordinates on S'^. We see that for sets B{0, 1 + R~^) C -R > 1 in original 
coordinates we can use the superharmonic potentials 

Mr.*.«)=mi„(l,i^) = ,„i„(l,|^) 

to see that the Brownian motion sent from x E Cl\D does not enter D a.s.. 
By writing the solution u in terms of spherical harmonics, we have 

oo n 

u{t,4>,9) = Y, E an,me-''-'Y:;,{4>,9) 

n=l m=— n 

where A„ = p~^^n{n + 1). Using the original coordinates (r, ^,6) of \ D 
we see 

oo n 

uir,4>,9) = J2 E an,m(r-l)^"y;:(<^,^). 

n=l m=—n 

Thus, if p is small enough, the solution goes to constant and its derivative 
to zero when r ^ 1"*" so fast that the solutions of the conductivity equation 
satisfy the equation in sense of distributions. 
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